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It is shown that the total energy of the static "field + particle" system, defined in the framework 
of classical, renormalized electrodynamics of particles and fields, depends in an unstable way upon 
the field boundary data. It is argued that this phenomenon may be also an origin of the unstable 
dynamical behaviour of the system (i.e. existence of "runaway solutions"). It is proved that a 
suitable polarization mechanism of the particle restores the stability, at least on the level of statics. 
Whether or not it restores also the full, dynamical stability of the theory is still an open question. 
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I. INTRODUCTION 

Classical electrodynamics in its present form is unable 
to describe interaction between charged particles, inter- 
mediated by electromagnetic field. Indeed, typical well 
posed problems of the theory are of the contradictory na- 
ture: either we solve partial differential equations for the 
field, with particle trajectories providing sources (given 
a priori !), or we solve ordinary differential equations 
for the trajectories of test particles, with fields providing 
forces (given a priori !). Combining these two procedures 
into a single theory leads to a contradiction: Lorentz 
force due to self-interaction is infinite in case of a point 
particle. 

There were many attempts to overcome these difficul- 
ties. One of them consists in using the Lorentz-Dirac 
equation (see Q,Q,0). Here, an effective force by 
which the retarded solution computed for a given par- 
ticle trajectory acts on that particle is postulated (the 
remaining field is finite and acts by the usual Lorentz 
force). Unfortunately, this approach leads to the so called 
runaway solutions which are unphysical. 

Various remedies have been proposed to cure such dis- 
ease, most of them just based on a fine tuning of bound- 
ary conditions. Unfortunately, such a tuning excludes 
physically interesting problems (i.e. circular motion) and 
the question arises if one can construct a theory which 
does not contain unphysical solutions at all. The au- 
thors believe that to achieve the above goal we should 
first gain a deeper understanding of foundations of the 
runaway behaviour. 

As a starting point of our analysis, we use an approach 
proposed by one of us in papers Q and 0. It consists 
in defining an "already renormalized" four-momentum of 
the physical system "particle(s) + fields" . Equations of 
motion are then derived as a consequence of the conserva- 
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tion law imposed on this object. We deeply believe that 
such an approach is a a correct realization of the Ein- 
stein's programme of "deriving equations of motion from 
field equations" and that a similar procedure should be 
applied to formulate the two-body-problem in General 
Relativity Theory. 

We show in the present paper, that the physical in- 
stability is inherently contained in the renormalization 
method used. More precisely: in the simplest renormal- 
ization scheme the amount of energy contained "in the 
interior of the particle" decreases when the external field 
surrounding the particle increases. This contradicts the 
stability of the model. As a remedy for such drawback 
we propose the polarizability of the particle. Numerical 
analysis of such an improved model shows validity of this 
proposal. 

In this paper we analyze the renormalized energy of 
the total "particle + field" system on the level of statics 
only, but the energetic instability discovered this way is 
obviously a reason for the runaway behaviour of the dy- 
namical system as well. Indeed, the price which must be 
paid for acceleration becomes negative. This observation 
is fundamental, in our opinion, to understand the physi- 
cal reasons for the runaway behaviour of the theory and 
in search for a remedy for this phenomenon. 

The paper is organized as follows. In Section ITU the 
renormalization procedure proposed by one of us in £| 
(see also |2|) is presented. Then a monopole particle in- 
side a fixed volume V is considered: we compute renor- 
malized energy of the system and vary it with respect 
to particle's position. Next, we assume that the particle 
assumes position corresponding to minimal value of the 
energy. In this way we obtain total energy of the system 
as a function of the field boundary data, imposed on dV . 
Finally, we analyze stability of the system under small 
changes of these data. Here, both the Dirichlet-type and 
the Neumann-type boundary problems are considered. 

The above general results are then applied to a case of a 
monopole particle closed in spherical box. We prove that 
such system is not stable. Then we consider a polarizable 
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particle. Here, the external field may generate a non- 
vanishing dipole momentum, which changes completely 
the energy balance. It turns out that for a Heaviside-like 
relation between the field and the dipole momentum it 
generates, the system is stable. This suggests a possible 
way to improve in the future our renormalization method 
and to avoid (maybe) also dynamical instabilities, mani- 
festing themselves in the runaway behaviour. 



II. THE RENORMALIZED FOUR-MOMENTUM 
VECTOR 

Full description of the renormalized electrodynamics 
was proposed in Q or 0. In the present Section we re- 
view briefly heuristic ideas that stand behind definition of 
the renormalized four-momentum of the dynamical "par- 
ticle + field" system. 

As a starting point of our considerations take an 
extended-particle model. This means that we consider a 
fully relativistic, gauge-invariant, interacting "matter + 
electromagnetism" field theory, which is possibly highly 
non-linear. A moving particle is described by a solution 
of the theory, such that the "non-linearity-region" (or 
the "strong-field-region" ) is concentrated in a tiny world 
tube W around a smooth, timelike trajectory Q. We as- 
sume that outside of this tube matter fields practically 
vanish and the electromagnetic field is sufficiently weak 
to be well described by the linear Maxwell theory. The 
four momentum of the total system "particle + field" is 
obtained by integration of a (conserved - due to Noether 
Theorem) total energy-momentum tensor %: 

V X = I , (1) 

over a spacelike hyperplane E. 

We assume, moreover, that this fundamental theory 
admits also a static, stable, soliton-like solution, which 
will be called a "particle at rest". Here, the strong-field 
region (interior of the particle) is assumed to be concen- 
trated around the straight line x =const. Let m denote 
the total energy (mass) of this solution. Due to rela- 
tivistic invariance, we have also a six parameter family 
of solutions obtained by acting with Poincare transfor- 
mations on the static solution. Each of these solutions 
may be called a "uniformly moving particle" . If the so- 
lution has been boosted to the four-velocity u\ and if 
T(u) denotes its energy-momentum tensor, then the to- 
tal four-momentum of this solution equals mu\ and we 
have: 

mux = J T\{u)da^ . (2) 
This leads to a trivial identity: 

V x = mu x + jf (I^-T^u))^ , (3) 



which becomes extremely useful in the following arrange- 
ment. We assume that the straight line which describes 
the "trajectory" of the second (uniformly moving) par- 
ticle is tangent to the approximate trajectory £ of the 
first (i.e. generic) particle at their intersection point with 
S. If K(R) C £ denotes the ball of radius R, which 
contains the strong field region of both solutions, but is 
small with respect to the characteristic distance of the 
external Maxwell fields, then we have: 

Vx = mu x + f {%\-T\{u))d(7 li + 

JY,-K(R) 

+ I (T\~T\(u))da,. (4) 
Jk(r) 

Our assumption about stability of the free particle (soli- 
ton solution) means that the last integral is negligible 
since inside the particle both solutions are very close to 
each other. But the first integral contains only contri- 
butions from external Maxwell fields accompanying both 
particles. This way we have proved that the following 
formula: 

Vx ~mux+ f (1\ - T\(u)) da^ , (5) 

JT,-K(R) 

containing only external Maxwell field surrounding the 
particle, provides a good approximation of the total four- 
momentum of the total "particle + field" system. 

The theory proposed in £| consists in mimicking the 
above formula in the point particle model. Hence, we 
consider solutions of Maxwell equations having a "delta- 
like" current corresponding to a point charge e traveling 
over a trajectory £. Such a solution is treated as an ide- 
alized description of external properties of the extended 
particle considered above. Denote by T the energy mo- 
mentum tensor of this solution. Of course, the uniformly 
moving particle, whose four-velocity equals u, is repre- 
sented in this picture by a boosted Coulomb field, and 
its energy- momentum tensor is denoted by T(u). If tra- 
jectories of both particles are again tangent with each 
other at their common point of intersection with S, then 
momentum (JSJ may be rewritten as: 

Vx ~mux+ [ {T\ - T\(u)) da^ , (6) 

JT,-K(R) 

because outside of the particle 1 reduces to T and T(u) 
reduces to T(u). The main observation done in Q| is that, 
due to cancellation of principal singularities of both T 
and T(u), the above integration may be extended to the 
entire S. More precisely, the following quantity: 

Vx := mux +P f {T\ - T\{u)) da, t , (7) 

is well defined ( "P" denotes the "principal value" of the 
integral). According to the discussion above, we inter- 
pret this quantity as the total four-momentum of the in- 
teracting system composed of the point particle and the 
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Maxwell field accompanying the particle. Consequently, 
we impose conservation of V as an additional condition. 
This implies equations of motion of the point particle as 
a good approximation of equations of motion of the true, 
extended particle. 

This approach has an obvious generalization to the sys- 
tem of many particles (see |^). Also polarizable parti- 
cles, carrying magnetic or electric moment (and - con- 
sequently - displaying stronger field singularity than the 
Coulomb field) may be treated this way (cf. §). Re- 
cently, the above approach was improved by replacing the 
reference Coulomb field in Q by the Born field, matching 
not only particle's velocity but also its acceleration. This 
way the principal-value-sign "P" may be omitted in the 
definition because the corresponding integral converges 
absolutely (cf. 0). 

In what follows, we are going to apply definition Q to 
static "particle + field" configurations only. 



we obtain the total energy contained in V: 

Uv = m- \ \ D sing 2 dv + \ ( D reg 2 dv + 

+ / D sing D reg dv. (12) 
Jv 

Given boundary conditions, we are going to minimize 
the above quantity with respect to the particle's position 
r*o € V. Assuming that the particle always tries to mini- 
mize the energy of the system, we can write both fo and 
the total "particle+field" energy as functions of the field 
boundary data. Stability of the energy with respect to 
the boundary data on dV will then be studied. Before 
we pass to the above programme, we must specify which 
kind of boundary conditions on dV have to be controlled. 

A. Neumann conditions 



III. ELECTROSTATICS OF A MONOPOLE 
PARTICLE 



Consider now electrostatic field D surrounding the par- 
ticle with charge e, situated at the point r*o Due to 
Maxwell equations, the Gauss law: 



VD = e8(r-r ) , 



(8) 



must be satisfied, where by S we denote Dirac delta dis- 
tribution (in contrast with conventional 5, denoting vari- 
ation of a function). It is, therefore, convenient to de- 
compose the field into its singular and regular parts: 



P D reg -\- D sing 5 



(9) 



where the singular part D S i n g is simply the Coulomb 
field: 



P 



sing 



e(r- f ) 
47r||r — r*o|| 3 



whereas the remaining field D reg := D — D s , mg is di- 
vergenceless: VP reg = 0. Moreover, static Maxwell 
equations imply the existence of the scalar potential <fi: 
D = — V0. Hence, we have: Ac/> reg = 0. 

According to J7J), the complete energy of this "particle 
+ field" system contained in the the entire £ equals: 



H = m 



[D 2 - D 2 smg ) dv 



(10) 



Varying the energy integral 1121) with respect to the 
particle's position we get: 



SHi 



/ {D reg ■ {5D reg + 8D sing ) + D sing SD reg } dv 
Jv 

(13) 



Dsing 8D reg dv. 



For Neumann conditions we put D = — V0 for both the 
regular and the singular parts of the field, outside of the 
variation <5. Integrating by parts and using V D reg = 
we get: 



SH V = 



(t>regS(VD sing )dv- [ {<j)5D^}dcJ. (14) 

Jav 

But the variation of iJSJ gives us: 

S(VD sing ) = S {eS (r- r )) = -ed k (5 (r - r )) 5x k , 

(15) 

where Sxq denotes a virtual displacement of the particle. 
Imposing Neumann conditions D^ldV — /, where / is 
a fixed function, we obtain: 5D 1 - = on dV . Hence, 
the surface integral vanishes. Inserting (|15fl into l|14|) we 
derive the following formula: 

SHv = -eDl cg (x^)Sxl (16) 

We conclude that the extremum of energy condition im- 
plies the following static equilibrium equation: 

Dl cg (x h Q ) = . (17) 



B. Dirichlet conditions 



We suppose that the particle is contained in a fixed vol- 
ume V 9 r*o . Subtracting from Ti the electrostatic energy 
contained outside of V: 



1 



D 2 dv 



(11) 



For Dirichlet case we put SD = — VScf) for both the reg- 
ular and the singular parts of the field and then integrate 
(f 131) by parts. We obtain: 



8H V = I {VD sing )5(t> reg dv- I {D 1 

JV JdV 



da. (18) 
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Imposing Dirichlet conditions <p\dV — /, where / is a 
fixed function, we obtain: 5(f> = on dV and, therefore, 
the surface integral vanishes again. To derive the equi- 
librium condition 117fl from the variational principle, we 
must perform the following Legendre transformation: 



(VD sing ) 54> reg dv = 51 (VD sins ) 4> reg dv+ 
v Jv 



(5VD sing ) 4>regdv . (19) 



Then we use (jSJ) and (|15p. This way we obtain: 
6 (H v - e<t> reg (r )) = Dl cs (x h )6x k . 



(20) 



Comparing l|ltj|) and (|2U[1 we observe that the equilibrium 
condition 117fl may either be obtained from the varia- 
tional principle 5 (Ti-v) — 0, when the Neumann bound- 
ary data are controlled, or from the variational princi- 
ple 6 (Ty) = 0, with Ty := Jiy — e4>reg{T$), when the 
Dirichlet boundary data are controlled. The quantity 
TLy is the total energy of the "particle + field" system, 
whereas Ty is an analog of the free energy in thermo- 
dynamics. We conclude that imposing Neumann condi- 
tion on the boundary corresponds to the adiabatic insula- 
tion of the system, whereas imposing Dirichlet condition 
means that we expose it to a kind of a "thermal bath" . 
Indeed, imposing e.g. condition <p\dV = we must cover 
the surface dV with a metal shell and ground it electri- 
cally. This means that we admit energy exchange of our 
system with the earth. Similarly as in thermodynamics, 
the free energy Ty, which we optimize, contains not only 
the system's energy Hy but also the term "— e<^ re9 (ro)" 
which we interpret as energy of the "boundary-condition- 
controlling device" . Of course, from the point of view of 
the particle, both conditions lead to the same equation: 
D re «(xg) = because our theory is local and the parti- 
cle interacts with its immediate neighbourhood only, no 
matter how the boundary data are controlled far away 
from the particle. 



IV. AN EXAMPLE - MONOPOLE PARTICLE 
IN A SPHERICAL BOX 



The problem consists in solving equation A<p = 
—eS(r — f*o), where r*o <= K(0,R). In the Neumann case 
we impose the following condition: 



D ■ ri\ _ = E • n 

I r—R 



AttR 2 



(21) 



where E is a fixed three dimensional vector. 

In the Dirichlet case we impose the following condition: 



D = -E-nR-\ 



(22) 



Because of the axial symmetry of the problem, we may 
restrict ourselves to the analysis of the energy functional 
at points fo which are parallel to E: fo\\E. With this 
simplification, we are able to find an explicit solution 

4> = (f>sing + 4>reg, where: 

1 e 



47T \r — ro | 



in both Dirichlet and Neumann cases (cf. Appendices 
and[nj. To write an explicit formula for <fi reg it is useful 
to introduce the following variable: 



ro 



\E\ 



<E\ro) 



which runs from —R to R. Under this convention we 
obtain: 



IJ reg 



R 



1 



4tt V Vi? 4 + r 2 r 2 - 2r rR 2 cos 6 R 



-R l - 



R 2 - r r cos 6 + \J R 1 + ro 2 r 2 - 2r rR 2 cos 9 



- Er + —\n(2R 2 ) , 
R 

in the Neumann case, whereas: 



(23) 



Iy reg 



R 



47r \R V-R 4 + r ' 2 r 2 - 2r rR 2 cosf 



- Er 



(24) 



in the Dirichlet case. 



A. Stability 



In this section we shall analyze stability of a charged, 
monopole particle closed in a spherical box with radius 
R: V = K(0, R) C R 3 . Simplicity of the model allows us 
to solve explicitly the static Maxwell equations (for both 
the Neumann and the Dirichlet cases) and to compute 
renormalized energy of the system. Then we will find 
the extremum of the energy function with respect to the 
particle's position and check that for the Neumann case 
we get the minimum and for the Dirichlet case - the max- 
imum of the energy. Assuming that the particle always 
minimizes the energy, we will express energy function in 
terms of the boundary data and show that the system is 
unstable under small changes of these data. 



In both cases, the renormalized energy can be com- 
puted explicitly. Denoting E := \\E\\ we obtain the fol- 
lowing result: 



R 



1 



2 V 4tt V R 2 - r 2 



-nR A E 2 - 2eEr ( 



1 - 



R 2 



in the Neumann case (cf. Appendix [BJ) and: 



R 



4tt R 2 



(25) 



(26) 
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in the Dirichlet case (cf. Appendix . Finally, we com- 
pute the electric "free energy" T = Ti — e(j> reg (r ) in the 
Dirichlet case: 



T = m 



1 e 2 R 4 , 9 e 2 2 
s T + 2e£;r + -it RE — 

2 \4nR 2 -rl 3 4tt i? 



(27) 

We see that the equilibrium condition in the Neumann 
case reads: 



I ^ = <^> eE 



'o 



4tt R{R? -rg) 
whereas in the Dirichlet case it reads: 

- 2 i?r „ 5 



= 



e_D r 



4tt (i? 2 - 



2\2 



or 



(28) 



(29) 



We express the energy in terms of the following, stan- 
dardized variables: 



To , r 4nR 2 



Denoting: 
we obtain: 



1 



1 -x 2 

1 2 



n = (H - m) — — 



-ln|l-x 2 | -2qx 
1 



l-x 2 



(30) 

(31) 

(32) 
(33) 



Observe that for q = both energies may be expanded 
as follows (cf. figure 2]): 



H'm = -l + 2x 2 + 0(x 4 ) 

n'v = -i - * 2 + o(x 4 ) . 



(34) 
(35) 



This implies that only in the Neumann case the equilib- 
rium point (x — 0) is also a minimum of the energy. In 
the Dirichlet case the energy has a local maximum at the 
equilibrium point. As may be easily seen, this happens 
also for any value of E, Hence, for the Dirichlet case the 
free energy T should be used, for which local extremum 
is also minimum. In what follows we shall use the local, 
physical energy and consequently, we restrict ourselves 
to the Neumann case only. 



FIG. 1: Graph of renormalized energy vs particle's position 
and q — for H! m and TC'd 



For small values of q this enables us to express equilib- 
rium position in terms of the boundary data: 



(37) 



The same result could be obtained from the following 
expansion: 



H'm{x, q) = -1 + \q 2 - 2qx + 2x 2 + 0(x 4 ), 



d x 7i'jsr{x, q) = =>• x 



1 



(. 



q 2 + 0(q 3 ). 



(38) 
(39) 
(40) 



Observe that for increasing values of q, the energy of the 
system decreases (cf. figure 0)! The system "particle + 
field" turns out to be unstable - even small fluctuations 
of the external field q can decrease its total energy. This 
means that the particle behaves like a perpetuum mobile, 
providing a source of energy at no costs. In our opin- 
ion this unphysical feature of the model, manifestly seen 
in its static behaviour, could possibly be a source of its 
dynamical instability, i.e. the existence of "runaway" so- 
lutions of Dirac equation. As a remedy, described in the 
sequel, we propose to equip the particle with an addi- 
tional mechanism which, via electric polarizability, will 
restore its static stability. 



B. Neumann conditions 



In terms of the standardized variables, the equilibrium 
condition (|28(l reads: 



x(2-x 2 ) 
(l-x 2 ) 2 



(36) 



V. POL ARIZ ABLE PARTICLE 

We assume that the particle may get a non-vanishing 
electric dipole moment due to interaction with the neigh- 
boring field. We prove in the sequel that, under a suit- 
able choice of the polarizability properties of the particle, 
the resulting "particle + field" system becomes statically 
stable. 
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FIG. 2: Graph of renormalized energy vs boundary field q for 
~H'M(x(q),q) 



For a polarized particle, formula l|12|l for the total en- 
ergy remains valid but the field singularity is now deeper 
than in J3J), namely: 

VD = VD sing = eS(r- f ) - P k d k S{f- f ) , (41) 

where p k is a dipole moment. We assume that p k has 
been generated by the surrounding electric field D ac- 
cording to some law p = p(D reg (ro)), describing the sen- 
sitivity of the particle. Moreover, we admit the depen- 
dence of the coefficient m in (|7|) (and, consequently, in 
$l'2\i ) upon polarization. It will be shown in the sequel 
that insisting in having m constant we are not able to 
make the model physically consistent. Moreover, it will 
be shown that the electric sensitivity is uniquely implied 
by the dependence m = m(p). 



A. Variational principle 

Variation of the renormalized energy (|12[) with respect 
to the particle's position contains now the non- vanishing 
term 5m. Similar calculations as for the scalar particle 
lead, in case of the Neumann boundary conditions, to 
formula: 



57iv = 5m 



5(VD aing )dv 



[ {<t>6D^} da , 

JdV 



(42) 



lav 

and, in case of the Dirichlet conditions, to: 



da 



5Hv = 5m + (VD 8ing )54> reg dv - I {D J 
Jv JdV 

5m + 5 / (\JD sing )4> reg dv - I 4> reg 5{\J D sing )dv+ 
Jv Jv 



I'dV 



da 



(43) 



According to (|41|l . the new version of formula (|15f) reads: 
5{VD smg ) = - (ed k S(r- r ) - p>d 3 d k d(r- r )) 5x k Q + 



(d k S(r-r a ))5p k . 



(44) 



Plugging (|44|l into l|42|l we see that the total energy vari- 
ation splits into the sum of two pieces: the work due to 
virtual displacement of the particle and the remaining 
work, due to variation of m and p: 



5H V = - (eD reg +p k d k D reg ) \ 5r + 



5m 



D reg\ f=fo 5p 



(45) 



The second part B is obviously nonlocal - both the mass 
m and the moment p depend upon the value of D reg (r$). 
This quantity must be obtained from the field equation: 



A0 r 



0, with boundary value depending upon the 



particle's position. The only way to save locality of the 
model is to force the term B to vanish identically by 
imposing the following constraint: 



5m = D reg (f )5p 



(46) 



Denoting by mo — m(0) the mass of the unpolarized 
particle and by f(p) the additional polarization energy: 



m(p) = m + f(p) 
formula (|46|) may be written as: 



D k 9 (ro) 



dm 

dp k 



(47) 



(48) 



We see that the polarization energy / must play role 
of the generating function for the polarizability relation, 
otherwise the model would not be local. Indeed, suppose 
that B does not vanish and the particle's equilibrium 
condition needs vanishing of the whole right hand side 
of l|4"5|l . To decide whether or not its actual position is 
acceptable as an equilibrium position, the particle must 
know not only the field in its immediate neighbourhood, 
but also the shape of V and the field boundary data on 
dV. Such a behaviour is physically non acceptable. 

Inverting the generating formula (|48|l . we may find the 
dependence p — p(D reg (r )), which is uniquely implied 
by the "equation of state" (|47|l . Hence, we have: 



5H X 



(eD reg +p~d k D re 



L _ 5f , 

\r—ro 



(49) 



and the equilibrium condition becomes a local equations: 



(eD reg +p k d k D reg ) |_. =r _ = 



(50) 



A similar procedure works in the Dirichlet case as well. 
Applying the state equation to (|43|) we obtain: 



5T V = (eD reg + p k d k D reg ) | _. = _, 5f 



(51) 
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where the "free energy" Ty is given as: 

Ty := TLy - I (y D sing)4>Teg ~ 2 / 
JV 



:= H v - e(f> reg (r ) + D reg \ ?=?a -p-2f . (52) 

Equilibrium condition SJ-y = reduces to the same, local 
equation JSOJ). 



VI. AN EXAMPLE - POL ARIZ ABLE PARTICLE 
IN A SPHERICAL BOX 

Let us come back to the simple model described in 
Section IIVI on page 21 For the polarizable particle we 
must solve the field equation: 



A0 = -eS{f- f ) +p- V(<J(r - f )) , 



(53) 



where r*o € K(0,R), with either Neumann l|21[l or Dirich- 
let condition (|22() . We want to compute renormalized to- 
tal energy of the "particle + field" system and to prove 
that for a suitable state equation l|47[) our model becomes 
stable. 

Splitting the solution cf) into two parts: 

= <prnon + ^dip ^ ^ 

where by (j> mon we denote the solution of the monopole 
problem, found earlier (cf. Section l|IV|) . page^J), we 
reduce the problem to equation: 



A 



p ■ V (S(f - f )) , 



(55) 



with homogeneous boundary conditions: D dlp •n\ r _ R = 
in the Neumann case and 4 >dlp \ r -n = in the Dirichlet 
case. Choosing the axis e z parallel to E and passing to 
spherical coordinates (r, 6, tp) we obtain for tq — (ro, 0, 0) 
and p — pe z + p x e x (see Appendix iDl on pagc lll(l : 



,dip 
reg 



where: 



.dip = 1 p ■ (r - r ) 

sing - 4?r \ft_jt Q \- - 



/dip 
reg 




(56) 



(57) 



R 3 (R 2 -rr cos0) 




r (R 4 + {r r) 2 - 2rr R 2 cost 
rR 3, sin 9 cos tp 



(R 4 + r 2 r 2 - 2r rR 2 cos 9) 5 

cos ip{R 2 cos 9 — r$r) cosflcost/A 
Rr sin 9^R 4 + r 2 r 2 - 2r rR 2 cos 9 Rr sin 9 J ' 

(58) 

As we already noticed in the monopole case, axial sym- 
metry of the problem implies that minimum of the energy 
is assumed at the point ro which is parallel to E. The 
same argument implies that we have p x = in this con- 
figuration. We are going to limit our analysis to such 
configurations only. 



A. Stability 

We compute the total, renormalized energy of the sys- 
tem as a sum of two parts: 



U = W 



U dip ^ 



(59) 



where T-[ mon denotes the energy of the monopole field 
obtained earlier f l|25|l . pageQJ, and T-i dl P denotes the re- 
maining part, containing energy of the dipole field and 
the interaction energy. The latter term is computed in 
Appendix lEl (page 1 1 3fl . The final result for the Neumann 
case, written in terms of standardized variables reads: 



l-T 2 

2 ( p x{2- x 2 ) 

3 [eRjl 



1 —-ln\l-x 2 \-2qx + ^q 2 



V 



1 



P 



) 2 e 2 R 2 (1 - x 2 ) 3 eR 



2+ 
(60) 



Now, stability of the system depends upon the polariz- 
ability of the particle, i.e. upon the choice of the "state 
function" / (cf. (|47(l on page EJ). At the moment we 
have no general criterion which would guarantee stabil- 
ity. However, it is easy to show that for: 



f(p) = -j\\p\\ 3 



D reg = ~c 2 \\p\\p, c>0, (61) 



our system is stable. Indeed, using l(2*3l and l(55|l we 
obtain the following equation for the value of the dipole 
moment p: 



c 2 p 2 sgn(p) = D 



reg 



= -Vi 



irnon 
reg 



l 



eg 



2p 



ex{2 - x 2 ) 

4tt \R 2 R 3 (l - x 2 ) 3 ~ R 2 (l~x 2 ) 2 



(62) 



Denoting Airec R = C and p = ^L, we get equation for 
p: 



-Cp 2 sgn{p) 



q 



2p 



x(2-x 2 ) 

(1-x 2 ) 3 (1-x 2 ) 2 



(63) 



For small x, we use Taylor expansion of the right hand 
side. Consequently, we have: 



-(7prsgn(p) w q — 2p — 2x — 6px 2 — 3x 3 



(64) 



For p > there are two solutions of this equation for 
small x and q: 



Pi 

P2 



c{ 1 + ^ B+ VT^U 



For p < there is only one solution for small x and q: 



(65) 
(66) 



1 



P3 



rC 



(67) 
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Inserting the above solutions into the energy function 
(|rJUj> we define for i = 1,2, 3: 

H'i(x, q) = TLh'{x 7 q, eRpi) . 

It turns out that Ji! i does not admit any minimum with 
respect to x (i.e. a stable "field + particle" configuration). 
For the remaining two cases we use Taylor expansion for 
small x: 



n ^- i+ r -sea 



i 



3 VI - gC 



C 2 C 



9 



2g 1 



2 /l 1 
VT = g77 U ~ C 2 



■>|l_i_ + ±_i^_ 
" 1 C 2 C 3 1 - gC* 



(68) 



1 9 

- 1+ 3* 



2g 1 



1 



2 f 1 _ J_ _ q 
y/TTqC V3 C 2 C 



3C 2 3 VI + gC 
2 

C 2 



_f__ 9, 2 



2 1 



g_ 

c 



l l 



3 1 



gC 
(69) 
(70) 



Minimizing both energies with respect to x we obtain: 
3C 2 /„ 2C(C 2 -9) 2 



xi(q) 



x 3 (q) 



32(C* 2 - 3) 
C 2 (2C 4 - 15C 2 



C 2 -3 



-Q 



{C 2 - 3) 2 
3C 2 /„ 2C(C 2 



32(C* 2 - 3) V C 2 - 3 

C 2 (2C 4 - 15C 2 + 45) 



{C 2 - 3) 2 

Plugging Xi(q) into the energy we get for small q: 



(71) 



(72) 



W'i(«)«-l-57« 



15 + 4C 2 2 

g - 



3C 2 6(3 -C 2 ) 
(18 + 42C 2 -7C 4 )C* , 



W's(g) w -l - 



12(3 -C 2 ) 2 

15 + 4C 2 2 



3C 2 6(3 -C 2 ) 
(18 + 42C 2 - 7C 4 )C 3 
12(3 -C 2 ) 2 9 ■ 



(73) 



(74) 



We see that for C G]0, V3[ the g 2 term is positive. This 
means that the system "particle + field" does not be- 
have any longer like a perpetuum mobile: to deform 
its original configuration, corresponding to q = 0, the 
boundary-condition controlling device must perform a 
positive work. Hence, the system is stable under small 
changes of q (see figure |3J). 




FIG. 3: Graph of Tt'(q) - renormalized energy vs boundary 
field for dipole particle, C = 1 



B. Conclusions 

We have shown that the polarizability of the parti- 
cle, described by a suitable "state function" / (e.g. by 
lED); ma Y be a good remedy for the static instability 
of the renormalized electrodynamics of point particles. 
Whether or not this will cure also the dynamical insta- 
bility, i.e. the existence of "runaway" solutions, is another 
question which we would like to study in the nearest fu- 
ture. 

At the moment the bifurcation phenomenon occurring 
near the ground state q = is worthwhile to study. Ob- 
serve that the point r*o = 0, corresponding to q = 
and described by the purely monopole field, is not sta- 
ble. This configuration corresponds to a local maximum 
of the energy and belongs to the unstable branch of sta- 
tionary points, described by the function H'a- 



APPENDICES 

APPENDIX A: NEUMANN SOLUTION FOR 
"PARTICLE + FIELD SYSTEM" 

We are looking for a solution of the Poisson equation 
A<p = — e5(f — r*o) with boundary condition (|21l) . where 
I! t^o II < R an d ro|| £7. Denote: 



t>sing + 0reg 



Er , 



(Al) 



where 



1 e 



smg — 4w \p-p \ > <-^Vreg 
1 



A<#> = and: 



D° eg -n\ r=R = n- 



\r-r 



r=R 



AttR 2 



(A2) 
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To find 4>®egi we use the following formula (cf. P-83): 



1 



-i = ^P„(cos^)-^, (A3) 



r ' 



\Jr 2 + r 2 — 2 rr cos 9 r 

(6 is the angle between f and E) valid for — r < ro < r, 
together with the following Ansatz: 



6° 

reg 



J2cnr n P n (cOs6). 



(A4) 



Write boundary condition as 
d 



e_d_ (\ 1 

r=R 4?r dr V \?- 7 "o| 



(A5) 



'///, 
v///« ;; -\\ / 

\ / /- 

////// A / -'—»->, \ \ < 

If/// / \ \ 1 

/////////// \ WW \ \ '. . 




FIG. 4: Directions of the field D - _B for it = 1, r = 0.5 



and substitute (|A3|> and ljA4() to (|A5|) . This way we get 
the solution given as a series: 



efl^/ 1^ (rr ) 



(i? 2 )' 



— P„(cos0) . (A6) 



Observe that (|A3|) gives, after rescaling, the first compo- 
nent of l|A6|l . The second one will be obtained from the 
following: 



Plugging R 2 instead of r and rrg instead of r in Lemma 
(jA.lj) yields: 



€eg = 

e ( R 11 

4^ V \fW+r 2 r 2 - 2r Q rR 2 cos 9 R R 
1 
R 



ln(2P 2 )- 



Tn 



R 2 — tqt cos 9 + \/ R 1 + ra 2 r 2 — 2rorR 2 cos ( 



(A12) 



Lemma A.l For ||ro|| > r we have 



V-^p ( 

n— 1 



cose/ 



1 



In 



1 d r ° i 
- 1 cos ( 

2 I r 



— 2 — cosf^ 
r / r 



Figure shows the directions of the field D — E = 
D S ing + -Dreg + Tk^ ^ • Observe that the field is tangent 
to the boundary of K(0, R). 



APPENDIX B: RENORMALIZED ENERGY FOR 
NEUMANN SOLUTIONS 



Proof: Substituting t for tq in (|A3|I : 



To compute integral (|13fl : 



oo n _ x \ 



\Tl=l 



1 



1 



\dt 



t\Jr 2 + t 2 - 2rt cos 9 tr 



n=l 
1 



1 r ™ 


n r n+1 




r 


0" 






7 



In 



1 I i r ° , 

— 1 cos ( 

2 \ r 



(A7) 
(A8) 
(A9) 



as + i + t 2 - 2rtcos9 + In 



(A10) 



— 2 — cos( 
r / r 



(All) 



V 



E* sinqE re _gdv , 



observe that: 

1 f 2 , 1 

-2 Ly D -" iV = 2 



(Bl) 



i>singD- L sing da , (B2) 



3V-c 



bregD^reg da . (B3) 



Integrals containing products of singular and regular 
fields are understood in the sense of distributions (cf. , 
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p. 748). Denoting k e := K(ro,e) we obtain: 



D sinqD r& ndv — lim / D s i n nD re ndv — 



V-k c 



lim — 

+0 2 .i v -k c 



^ {^sing^reg H~ 4^reg^ sing) dv 



= — lim — 

e^O 2 jQ V -dk e 



^singD reg 4*regD sing) d(J . 



(B4) 



Knowing <j) we can compute TL: 

<j)\ r=R = -ERcos9+ 
e ( 2 



(Bll) 



1 1 

4vr V y/R 2 + r 2 - 2r Rcos9 R + R 



\n(2R)- 



In 

R 



R - r cos 9 + y/R 2 + r 2 - 2r Rcos9 ) , 

(B12) 



Hence, for V = K R := K(0, R) we have: 



U=m- - I (pD 1 - da+ 



lim 



1 

dk. 



5 2 yyreg-^stng 



J s%ng 



Dreg) da. (B5) 



The formula is true for both the monopole and the dipole 
singularity of D S i ng . Here, we consider the monopole 
(Coulomb) singularity. In this case the function 4> S i ng 
multiplied by e 2 (coming from the surface measure da) 
vanishes for e — > 0. Hence, we have: 



U = m-- [ <\>D^ da + lim - / (f) reg D± da . 
2 JdK R e ->° 2 J dkt 

(B6) 

To compute the integral over dk e , we use spherical coor- 
dinates (e,P,cp) centered at ?q. Parameters r and cos# 
present in 4> reg may be expressed as follows: 



r 2 = Tq + e 2 — 2ero cos r cos 6* = tq — e cos j3 



lim r 2 = rn , lim r cos = rn 

e->0 U e->0 



(B7) 
(B8) 



Then: 



-e£> + — 



1? 



i-im 



Note that: 



4tt \R 2 -r 2 i? i? 



zdcr = AttR, 



1 - — 

i? 2 



A' R V-R 2 + t"o 2 - 2r i?cos 

E cos 9da = 0, 
if 8 



(B13) 

(B14) 
(B15) 



/ In R — rn cos 9+^R 2 +r 2 - 2r Rcos9 

J Kr 



da = 



= 4irR 2 In 2R, 
where we used two integrals 2.736 from 0. Then: 



4ttR 2 



Kr 



AttR 



Moreover: 



(B16) 



(B17) 



E Z R I cos 2 (Mr = --kR a E\ 

Kr 3 

E cos 9 ,4 . . 

drr = -nEr , (B18 
Kr VR' 2 + r 2 -2r Rcos9 3 

/ In R - r cos 9 + \/R 2 + r 2 - 2r i?cos( 

JKr 



'Kr 

4 

x E cos 9da = — -7rrni?, 

O 



(B19) 



1 



lim - / cj> reg D si da 
e27r r 1 



lim 

e^O 2 47T J e 



2</ , re g ('"o,e,/3) sin /3e 2 d/3 



-r^>reg(ro,r = r ,9 = 0) / sin/3d/3 = -e 



1 



2 ^ re 9 I r=r 



(B9) 



where we used four integrals 2.736 from |l2j. Then: 

E [ <j) cos 9da = --itR 3 E 2 + — ( -7rr + -nr 
Jkr 3 4tt V3 3 

= --irR 3 E 2 + eEr . (B20) 
3 

The final result is the sum of JBT3j|, l|B16)l and (|B20|) 

with coefficient i; 



Consequently: 



TL = m 



dK R 



4>D da H — e<t> rPn L _ . 

r 2 ^ ' \r—ro 



(BIO) 



TL =m 



R 



2 V 4tt V R 2 - 



-nR 3 E 2 - 2eEr 



R 



In 



R 2 



(B21) 
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APPENDIX C: DIRICHLET SOLUTION AND 
THE CORRESPONDING ENERGY 



APPENDIX D: DIPOLE PARTICLE IN A 
SPHERICAL BOX 



To find a solution of the Poisson equation A<j) = 
— eS(r — fo) with boundary conditions (|22[l . where ||fo|| < 
R and r \\E, we denote: (p = 



47T f-f l ' " 



Wj + 0re 3 - -Er, where 



and: 



reg 

47r V f — r*o 



r=,R 



4?ri? 



(CI) 



Again, we use Ansatz (|A4|) as we did in Appendix ^ 
page|Hl and expand also boundary conditions: 



6° 

reg 



en 

47r V r 



\r-r 



(C2) 



r=R 



in series of Legendre polynomials. After substitution 
(lAl)l and (53} to fC^jl we obtain: 



6° 

reg 



eR x ■> 
1^ 2^ 



(n"o) 



4- ^ 

After rescaling (|A3ft we get: 



-^■PnicosO). 



reg 



l 



(C3) 



(C4) 



4tt \R V-R 4 + r 2 r 2 - 2r rR 2 cos , 

Singular part of the electric field has the Coulomb singu- 
larity at r*o- Hence, formula (|B10|) is valid. However, we 
have: 



<f>\r=R 
e 4>reg\ f= p o 

This implies 
2ttR 2 



R 2 



4irR (B? + ro 2 _ 2r Rcos6p 



-ER cos 6 



e 1 
47 



-e£r + |- f - 



1? 



4?r V i? i? 2 - r 2 



.E cos 6, 
(C5) 
(C6) 

(C7) 



(i? 2 



4) sin8 d8 



(4W?) 2 (i?2 



r 2 — 2r o i?cos0)'2 



47ri?' 
(C8) 



-2tt.R 2 



-2?ri? 2 



Si? cos ( 



{R 2 - r 2 ) sin 6 d6 



4irR(Ri +r 2 -2roRcos9) 

-eEr 0l 

E 2 R cos 2 9 sin 6 dO = - -irR 3 E 2 , 



cos 6 smOdB = 0. 



Consequently, we obtain: 
H = m+ I (^ttR 3 E 2 



2 V 3 



R 



Air R? - 



or, in standardized variables (|30fl . 



l-x 2 



(C9) 
(CIO) 

(Cll) 
(C12) 
(C13) 



We must solve equation A0 dlp = p • V (£(f — fo)) with 



boundary conditions D dip 



.dip 
sing 



\r=R 



0. Denoting 



reg 



where 



,di P _ . 1 p- (r - f ) 

sing 



47r |r — rop 



(Dl) 



we get Laplace equation A0* p = with boundary con- 
dition: 



Z3 dlp • f I 



ra • — V 

47T 



f (r - f ) 
If - f | 3 



(D2) 



r=_R 



For any pair of vectors f i p we choose coordinates in 
which fo is parallel to the z-axis e z and polarization vec- 
tor assumes the form p — pe z + p x e x . The final solu- 
tion will be the sum of two harmonic functions fulfilling 
boundary condition (|D2Jl . calculated separately for pe z 
and p x e x . 



Observe that, for 



imon 
J reg 



Vo,f) being a solution of 



Laplace equation, also the function | • Vf 4>W a" ^ s nar ~ 



monic. Moreover, if <^>™°™ fulfills conditions 



5J condi- 



tion from page El): 



|:C 9 Wo) 



r=JJ 



47r dr \r 



\r-r \ 



r=R 



(D3) 

then, after differentiation with respect to fo we obtain: 



o \ J r Vreq 

or \ e a 



— — V 1 

AirdrV' ?0 |f-fo| 



(D4) 



r=R 



Hence, the function | • V^ o 0"eg™ satisfies boundary con- 
ditions (ID2I). We conclude that: 



1 

47re 



(p'V ?0 )c: 



(D5) 



(cf. pl|. p. 14). Applying (|D5(1 for p = pe z + p x e x al- 
lows us to solve the problem separately for p parallel and 
orthogonal to fo. 
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1. Solution for p\\ r*o 

To obtain the parallel part we differentiate monopole 
solution f (|A12|) . Appendix |A"|) along the e z -axis: 

p d 



idip £^ " rrnon 

reg ~ e dr 9re9 



p d 



R 



4tt dr V VR 4 + r 2 r 2 - 2r rR 2 cos 9 R R 
1 



4 + 4 ln ( 2 ^ 2 ) 



-In 
R 

P 



R 2 - r r cos 9 + V ' R A + r 2 r 2 - 2r rR 2 cos ( 
R(r Q r 2 -rR 2 cos 9) 



4tt V (i? 4 + r 2 r 2 - 2r rR 2 cos 0) I 
1 

cost/ 

-r cos 9 (V-R 4 + r 2 r 2 - 2r rR 2 cos 6> + i? 2 ) + r 2 r 



i? 2 - rr cos 6* + \/i? 4 + r 2 r 2 - 2r rR 2 cos 6* 



But: 



(i? 2 - rr cos + + ^o 2 r 2 - 2r ri? 2 cos 6») 



(D6) 



(D7) 



x yR 2 - rr cos - \J R i + r 2 r 2 - 2r rR 2 cos ( 
= -(r r) 2 sin 2 9, (D8) 
(-rcosf? (V-R 4 + r 2 r 2 - 2r ri? 2 cos 9 + i? 2 ) + r 2 r 
x ( R 2 — rr cos 9 — \/i? 4 + r 2 r 2 — 2r rR 2 cos ( 



-r r 



sin 2 6»\/-R 4 + n> 2 r 2 - 2r ri? 2 cos 9+ 



r cos 9(—r rR 2 cos 0) + R 2 r 2 ro 



= r 2 ro sin 2 



ro 2 r 2 — 2tqtR 2 cos( 



(D9) 



So: 



= JP_ ( _ R(r r 2 -rR 2 cos 9) 
reg 4vr V (i? 4 + r 2 r 2 - 2r ri? 2 cos 9) I 



r 2 r sin 2 6> (i? 2 - V# 4 + f 2 r 2 - 2r rR 2 cost 
R(rr ) 2 sin 2 9 ^R 4 + r 2 r 2 - 2r Q rR 2 cos 9 
if pR 3 (R 2 - rr cos t 
4n \r (i? 4 + (r r) 2 - 2rr R 2 cos 




(D10) 



Figure [5] shows the directions of the field D dlp . Observe 
that the field is tangent to the boundary of K(0,R). 



2. Solutions for p_lr~o 



For p = p x e x we get: 



kdip _ . 
J reg — fx 



v raVreg 



Px d ^ 

e feo 



mon 
reg ■ 



(Dll) 



////// 

////// / 

/// i w-— -/ 


■-/ / ( I I \ \\\ 

W^^y/ i \ 


l \ \^~— ^/ / 1 
iVw^— - -// ! 

\ — ■ s/ 1 I 

\ v^-~- --v / 
i uW/ / 

\\\\ / 


\ / 

\ N~=^ / / / ( 



FIG. 5: Directions of the field D%* g + D*f g , R=l, r = 0.5, 
p = 1 



The easiest way to calculate this derivative is to use 
spherical coordinates ro = (ro, #0; Then: 

x = r Q sin 9 sin V3 &, yo = r sin 6» cos tp , 

zq = rocos^o , (D12) 



and: 



- — = sin 9 cos (po — - 

sin <po d 



a 1 9 

COS 9q COS (fi — — — - 

r d6» 



r sin 6»o dip 



(D13) 



But for rb||e z this procedure is singular because sin#o = 
0. To overcome this difficulty we first calculate the result 
for r*o !/f e z and then pass to the limit 6q — * and ipo —> 

0. For this purpose we must be able to differentiate the 
function cos 7, where 7 is the angle between r and r*o, 

1. e.: 

r ■ 7% = rry cos 7 , (D14) 

or, equivalcntly: 

cos 7 = cos 9 cos 9q + sin 9 sin 0q cos((^ — cpo) ■ (D15) 
Hence, (ID13|) gives us: 
i) 



dx 



cos 7 



— cos 9 cos ifo (— cos 9 sin #0 + sin 9 cos 9q cos((p — (po)) 
ro 

sinipo sin#sin(<^ — ipo) ) — sin^cos^. 

To J fo-^o r 

(D16) 



This method allows us to calculate effectively the deriva- 
tive of the monopole field from Appendix 1X1 (u. |SJ along 



13 



Xq. The final result reads: 

rR 3 sin 9 cos <p 



dip _ Px_l 

re9 ~~ 4iry (R4 + r 2 r 2 - 2r rR 2 cos 0)1 
cos (fi(R 2 cos 9 — r$r) 



cos 9 cos ip 



Rr sin 9y/R 4 + r ' 2 r' 2 - 2r rR 2 cos 9 Rr sin 9 I 

(D17) 

We stress that the above function is regular at 9 = due 
to cancellations between the second and the third term. 



APPENDIX E: RENORMALIZED ENERGY OF A 
DIPOLE PARTICLE 

To calculate J{ dl P we use results of Appendix [Bj It 
turns out that in formula (|ij5|) . only the following non- 
vanishing terms were not taken into account in Ti mon : 



H di P = __ 



- [ <p di PD^da 

2 JdK R 



mon 
si7ig 



b dvp D ± 
sing reg 



da 



(El) 



where: 



idip 
reg 



1 

4-7T 





1 


r sing 


47T 


irnon 


e 


rreg 


4tt 


-Im 

R 


R 2 



p R 3 (R 2 — rr Q cos 6*) p 
s r (R 4 + (r r) 2 - 2rr R 2 cos 6) * r o R 

p(r cos 9 — ro) 



(E2) 
(E3) 



I? 



1 



4tt V V-R 4 + r 2 r 2 - 2r ri? 2 cos 9 R 



R 2 - r Q r cos 9 + V 'i? 4 + r 2 r 2 - 2r ri? 2 cos ( 



- £r cos 6>+ -ln(2i? 2 ), 
ii 

= 0reg + 0smg ■ 



Moreover, we have: 



1 e 

— — + #cos( 

47T i? 2 

-— l(h mon - 

<j £ \Vreg 



(E4) 
(E5) 

(E6) 
(E7) 



To compute the integral over dK^ we note that: 

f,dip _ 

i? 2 -r 2 



-dip , .dip 



V 



reg r stng ^ , {R2 ^ ^ Q) g 



R 



(E8) 



whereas D 1 - is expressed by (|E6|) . Moreover: 
— 4t2tt / (j) dl P sin9 d9 = 0. 



(E9) 



So: 



—2ttE 
2 



b dip cos 9 sin 9 d9 = — pE. 

2 



To find the limit: 

lim — 

+0 2 ldk c 



b di P D ± 

sing reg 



(E10) 

we analyze behaviour of fields (|E2f) - (|ET|) for e — * 0. All 

these terms have at most the e _2 -singularity. Therefore, 
they are continuous and bounded when multiplied by e 2 . 
Thus, we can interchange the limit and the integration 
operations. 

We follow our procedure described in Appendix [B] 
page Using l|B7l) and l|B8|l we obtain in terms of the 
standardized variable x = ^: 



lim(e 2 ^J =-£ cos 0, 
d 



(Ell) 



lim 



imon 

~de 9 



r R 



ro 



4vr \(R 2 -r 2 ) 2 [R 2 -r 2 ) 2 R 
e 1 x(2-x 2 ) 
i^R 2 (1-x 2 ) 



E ) cos/3 



cos (3 — E cos j3 
-cos [3D 1 ; 



yinon 
reg 



lim 



d 



idip 



1 



2pR 3 



Se rre V 4n(R 2 -r 2 ) 3 
1 2p 



cos/3 



4tt R 3 (l ~ 
lim (e 2 D" lon ^ 



r 2\3 



cos (3 = - cos P DfP g 



si7ig ) 



e 

Itt' 



4tt i?(i? 2 



2\2 



p x(2 — a; 2 ) 
4^i? 2 (l -x 2 ) 2 

cos 2 /3 sin /3 d(3 ■■ 



> a 
Then: 

7 



r=r 

(E12) 



(E13) 
(E14) 

(E15) 
(E16) 



DTngKZ) da 



= I I l^Z. I n mon 

2 I 3 4^ 1 reg 



r>dip 
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Using q = ^-E and (jSH we obtain: 
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